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Abstract 

A special-purpose, semi-analytical solution method for detennining the stress and 
deformation fields in a thin laminated-composite cylindrical shell with an elliptical cutout is 
presented. The analysis includes the effects of cutout size, shape, and orientation; non-unifonn 
wall thickness; oval-cross-section eccentricity; and loading conditions. The loading conditions 
include unifonn tension, uniform torsion, and pure bending. The analysis approach is based on 
the principle of stationary potential energy and uses Lagrange multipliers to relax the kinematic 
admissibility requirements on the displacement representations through the use of idealized 
elastic edge restraints. Specifying appropriate stiffness values for the elastic extensional and 
rotational edge restraints (springs) allows the imposition of the kinematic boundary conditions in 
an indirect manner, which enables the use of a broader set of functions for representing the 
displacement fields. Selected results of parametric studies are presented for several geometric 
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parameters that demonstrate that analysis approach is a powerful means for developing design 
criteria for laminated-composite shells. 

Introduction 

Cutouts in cylindrical shell-type components are unavoidable in the construction of 
aerospace structures. This fact is significant because the structural failure of these components 
usually begins near the cutout because of high stress concentrations that initiate the fonnation of 
cracks. Hence, a cutout can trigger a local failure at a load level lower than the global failure 
load of a corresponding shell without a cutout. As a result, preliminary-design sizing of a 
cylindrical shell with a cutout is often based on the magnitude of the stress concentrations near 
the cutout. Therefore, an accurate assessment of the stress concentrations in a given shell 
subjected to various types of loading and support conditions is essential to the development of 
safe and reliable designs. Moreover, validated special-purpose analysis tools that enable rapid 
parametric studies would be very valuable to structural designers and for the development of new 
design criteria and design concepts. 

Several analytical, numerical and experimental studies have been conducted during the past 
sixty years to determine stress distributions in cylindrical shells with a cutout and subjected to 
various types of loadings; such as, axial tension and compression, torsion, and internal and 
external pressure. Pioneering analytical work was conducted by Lurie “ to investigate the effect 
of axial tension and internal pressure, and shell curvature, on the stress concentrations around a 
circular cutout in the 1940s. Many years later, analytical studies were presented by 
Lekkerkerker, Van Dyke, Ashmarin, Murthy et al., Guz et al., and Van Tooren et al. that 
further investigated the effects of various factors on the stress concentrations around a cutout in a 
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cylindrical shell. Similarly, experimental investigations have been conducted by Tennyson, 9 
Starnes, 10 Pierce and Chou, 11 Bull, 12 and Zirka and Chernopiskii, 13 and numerical studies have 
been conducted by Liang et ah, 14 and Shnerenko and Godzula. 15 In 1964 and 1972, respectively, 
Hicks 16 and Ebner and Jung 17 summarized the results obtained from several of these previous 
studies and provided extensive lists of references related to this problem. Most of these previous 
studies are for isotropic cylindrical shells with a circular cutout. Only a few of these studies, such 
as those presented by Pierce and Chou 11 and by Murthy et ah, 6 address the effects of cutout shape 
(elliptical cutouts) on the stress concentrations. 

Mitigation of high stress concentrations by tailoring shell- wall thickness, material orthotropy 
and anisotropy, and cutout reinforcement are also important considerations in the design of 
aerospace structures made of lightweight composite materials. Likewise, the potential for using 
shells of non-circular cross section are relevant to fuselage-like structures. However, only a few 
studies have considered these effects. For example, the influence of wall-thickness variation on 
the stresses in axially loaded composite cylindrical shells, without a cutout, has been investigated 
by Li et al. Although numerous analyses exist in the literature on the analysis of shells with 
circular cross sections, only a few include non-circular cross sections. Sheinman and Firer 19 
provided an analytical investigation of stresses in laminated cylindrical shells with arbitrary non- 
circular cross sections. More recently, Hyer and Wolford” " studied the effect of non-circular 
cross sections on damage initiation and progressive failure in composite cylinders by employing 
the finite element method. 

The objective of the present study is to present a special-purpose analysis for a laminated- 
composite cylindrical shell with an elliptical cutout that can be used to rapidly, and 
parametrically, investigate the effects of shell curvature; cutout size, shape, and orientation; and 
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ply lay-up on stress-resultant concentrations near the cutout. The analysis is applicable to thin- 
walled cylindrical shells with non-uniform wall thickness, a non-circular (e.g., oval) cross 
section, and subjected to tension, torsion, and bending loads as illustrated in Fig. 1. 

To accomplish this objective, an overview of the analysis is presented first. Next, the 
boundary value problem is defined along with the kinematics and stress-strain relations used in 
the analysis. Then, the derivation of the equations governing the response and numerical 
procedure are described. Finally, selected numerical results for oval and circular cylindrical 
shells with either circular or elliptical cutouts and subjected to either tension, torsion, or pure- 
bending loads are presented. 

Analysis Overview 

The analytical approach used herein permits the determination of the pointwise variation of 
displacement and stress components. It is based on the principle of stationary potential energy, 
but utilizes local and global functions that are not required to satisfy the kinematic boundary 
conditions directly. Thus, the choice of local and global functions is not limited by a particular 
type of kinematic boundary condition. The kinematic boundary conditions are imposed by 
employing the Lagrange multiplier method. Both local and global functions are used, in contrast 
to the traditional approach, to enhance the robustness of the analysis method. In particular, the 
local functions are used to capture rapidly varying stress and strain gradients and local 
deformations near a cutout. Toward that goal, Laurent series are used for the local functions and 
are expressed in tenns of the mapping functions introduced by Lekhnitskii. ~ Fourier series are 
used for the global functions and are used to capture the overall deformation and stress fields. 
The kinematic admissibility requirements on the local and global functions are relaxed by 
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defining that the edges of the shell are supported by extensional and rotational springs. Zero- 
valued displacement and rotation kinematic boundary conditions are enforced in an indirect 
manner by specifying values for the spring stiffnesses that are large compared to the 
corresponding shell stiffnesses. This approach effectively yields a prescribed kinematic 
boundary condition in the limit as the relative stiffness of the spring becomes much greater than 
the corresponding shell stiffness. Similarly, values for the spring stiffnesses can be selected that 
correspond to a given unifonn elastic restraint along an edge, similar to that provided by an end- 
ring. This capability is important, and useful, because in some test fixtures or actual structures 
the edge supports may not be stiff enough to simulate a fully clamped boundary condition or 
flexible enough to simulate a simply supported boundary condition. 

As suggested by Li et al. ls and Sheinman and Firer, 19 nonuniform wall-thickness variations 
of a shell, which lead to non-uniform laminate stiffnesses, are represented by using trigonometric 
series. Specifically, nonuniform shell-wall thickness is represented in the present study by 
perturbing the ply thicknesses with a function that is periodic in either the longitudinal or the 
circumferential direction. The variation in wall thickness is accounted for by adjusting the 
lamina properties, resulting in nonhomogeneous in-plane and bending stiffness matrices. The 
nonuni form shell curvature associated with a noncircular cross section is represented by using 
trigonometric series for the coordinates of an oval-cross-section shell reference surface. The 
aspect ratio, or out-of-roundness, of the cross-section is represented in the analysis by using an 
eccentricity parameter. 

In the derivation of the equations governing the response, the total potential energy consists 
of the elastic strain energy of the shell, the elastic edge restraints and the potential energy of the 
applied loads. The conditions that may arise from the choice of displacement approximations 
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without any kinematic restrictions are treated as constraint equations, and the potential energy 
arising from constraint reactions is invoked into the total potential energy through the use of 
Lagrange multipliers. The equations governing the shell response are obtained by enforcing the 
requirement that the first variation of the total potential energy vanish. The evaluation of the 
area integrals appearing in the potential energy are achieved numerically by using a basic 
quadrature method in conjunction with standard triangulation of the entire domain described by 
Shewchuk. Solution to the equations governing the response are obtained by using a standared 
Gaussian elimination procedure, which yields the generalized displacement coefficients and, 
thus, the stress and strain fields. The accuracy of the analysis depends on the number of terms 
used for the functional representation of the displacement fields. As the number of terms 
increases, the results converge to the exact solution. 

Representation of Shell Geometry 

The geometry of a thin-walled, noncircular, cylindrical shell of length L and with an elliptical 
cutout located at the shell mid-length is shown in Fig. 1. The origin of the global Cartesian 
coordinate system, (x, y, z ) is located at an end point of the longitudinal axis of the shell. As 
shown in Fig. 1, the x-axis coincides with the longitudinal axis of the shell. The y and z 
coordinates span the cross-sectional plane. A curvilinear coordinate system is also attached to the 
mid-surface of the cylindrical shell. The coordinates of points in the longitudinal, circumferential 
(tangential), and nonnal-to-the-surface (transverse) directions of the shell are denoted by (si, si, 
sj), and the corresponding unit base vectors are {ei, e 2 , e 3 } . 

Following Romano and Kempner, the non-circular cross-section of the cylindrical shell is 
defined as an oval with the coordinates, y and z expressed as 
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y = R 0 X a m(^) si 
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(la) 
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z = X & m(£)< 
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I cos — - 

m=l,3,5,7,... R o 


(lb) 


where E, represents the eccentricity of the oval cross section and R 0 is the equivalent radius of a 
circular cylindrical shell that has the same circumference as that of the oval cylindrical shell. 
The circumferential coordinate, s, varies between 0 and 27iR lt . The derivation of Eqs. (1) along 

25 

with the explicit forms of the coefficients aJJE) and bjc) are given by Madenci and Barut. 

As derived by Romano and Kempner, 23 the coordinates y and z in Eqs. (la) and (lb) can be 
related to the radius of curvature of an oval-cross-section cylindrical shell, R(s 2 ,£) , by 


R(s 2 ,& = 


R„ 


l + £cos(2 s 2 /R 0 ) 


( 2 ) 


Therefore, E, = 0 implies no eccentricity and corresponds to a circular cross section with 
radius R 0 . For positive values of the eccentricity parameter, c , the z -coordinate becomes the 
major axis and the y -coordinate becomes the minor axis. For negative values off;, the major 
and minor axes switch to the y- and z-axes , respectively. The range of values of the 
eccentricity parameter, E , , is bounded by-1 <£<1. 

As shown in Fig. 1, the cylindrical shell contains a cutout. The shape of this cutout is 
defined such that if the shell is cut along a generator and flattened into a plane, the cutout 
becomes an ellipse with major and minor axes denoted by a andb , respectively. For simplicity 
and convenience, the cutout is referred to herein as an "elliptical" cutout. Because the domain of 
the analysis shown in Fig. 2 corresponds to a similar flat region, a subsequent mapping of the 
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ellipse to a unit circle is possible, which enables the use of Laurent series expansions for the 
local functions. Note that the special case of a "circular" cutout is given by a = b . 

In the flat analysis domain, the minor and major axes of the ellipse are aligned with a local 
coordinate system, (x,,x 2 ) , whose origin is located at the center of the cutout and coincides with 

the origin of the parameter grid, given by constant values of si and S 2 , that forms the curvilinear 
coordinates (si, S 2 ) on the cylindrical shell mid-surface. The orientation of the elliptical cutout is 
arbitrary with respect to the longitudinal shell axis. Hence, the orientation of the local x, - axis 

(major axis) of the cutout and the longitudinal .s, -axis of the cylindrical shell is denoted by the 
angle, y/ . The elliptical coordinates, a and /? , representing a family of confocal ellipses and 
hyperbolas, respectively, are utilized in order to obtain the stress-resultant distribution in the 
direction tangent to the cutout boundary. The coordinate a is equal to a 0 = tanh '(b/a) on the 
particular ellipse that corresponds to the elliptical cutout. The other coordinate, J3 , varying 
fromO to 2k , is known as the eccentric angle and is related to the (x,,x 2 ) coordinate system by 
Xj = a cos /? and x, = b sin (3 . The eccentric angle /? is similar to the angle used for polar 
coordinates. 

The symmetrically laminated cylindrical shells considered herein are made of K specially 
orthotropic layers, and each layer has an orientation angle, 6 k , that is defined with respect to the 

5, -axis. Each layer also has elastic moduli E L and E r , shear modulus, G /r and Poisson’s 
ratio v i r , where the subscripts L and T represent the longitudinal (fiber) and transverse 


principal material directions, respectively. 



As for the shell thickness variation, the non-uniform wall thickness of the shell is denoted 
by h(s l ,s 2 ) , and its variation is included by assuming that the thickness of each ply, t k , varies as 
a function of the curvilinear coordinates in the form 


t k Oj , s 2 ) = f*o I 1 - £ l CoS I 2 T' \~ £ 2 CoS 


where t k0 denotes the nominal thickness of the k th layer in the laminate, and the parameters 


( m l , m 2 ) and (e 1 ,e 2 ) respectively, denote the wave numbers and the amplitudes of the periodic 


thickness variation in the longitudinal and circumferential directions. While the wall thickness of 
the shell is allowed to vary across the shell surface, the aspect ratio of the plies through the 
thickness is maintained, thus making the thickness variation of each ply to remain conformable 
to each other throughout the shell surface. A periodic thickness variation in the longitudinal 
direction is obtained by setting s x ^ 0 and s 2 = 0 , and in the circumferential direction by s x = 0 


and£' 2 ^ 0 . A shell with uniform thickness, t k = t k0 , is obtained by setting a\ = 0 ande 2 = 0 . 


( 3 ) 


Boundary Conditions and External Loads 

To facilitate a general imposition of prescribed boundary tractions, displacements, or 
rotations; the external as well as the internal edge boundary T of the shell is decomposed into 


r=r +r +r 

1 1 UW 1 (2 A 1 (3) 


( 4 ) 


As shown in Fig. 1, r (1) and f (2| denote the external edge boundary of the cylindrical shell and 


r (3) represents the traction- free internal edge boundary around the cutout. The unit vector 


nonnal to an edge is represented byn . Throughout this paper, a variable with the superscript 
is treated as a known quantity, arising from the externally applied loads or from prescribed 
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displacements and rotations. Also, the subscripts n, s, and t denote the directions normal, 
tangent, and transverse (through-the -thickness) to the boundary, respectively. The details of how 
prescribed edge loads and displacements are imposed in the analysis are presented subsequently. 


Prescribed edge loads 

External loads are applied to a shell by specifying values for the positive-valued stress 
resultants shown in Fig. 1. More precisely, the membrane loads applied to the f th boundary 
segment, T (<) , are given by 

N u =tl (5a) 

N n = t] (5b) 


where Nn and Nn are the axial and shear stress resultants, respectively, defined in the 
cylindrical coordinate system. Likewise, shell-wall bending loads that are applied to the t h 
boundary segment are given by 

M n = —m n (6a) 

M iu + 2M i2,2 = 6 - 2m* 2 (6b) 

where M n and M n are the pure -bending and twisting stress resultants, respectively, defined in 

the cylindrical coordinate system. Moreover, the left-hand side of Eq. (6b) is the Kirchhoff shear 
stress resultant of classical shell theory. 

As a matter of convenience, the analysis is formulated to also pennit the specification of 
concentrated forces and moments that are transmitted to the ends of the shell as if through a rigid 
end-ring, as shown in Fig. 3. Presently, the concentrated force P* and the concentrated axial 
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torque P* are included in the analysis. The force P* is simulated in the analysis by specifying a 
uniform distribution of the axial displacement, with the unknown magnitude A n , such that 


1 N, ,dT = p; 


(7a) 


Likewise, the torque P* is simulated by specifying a uniform distribution of the tangential 
displacement, with the unknown magnitude A , such that 


N u dT = P s 


(7b) 


The analytical process that is used to ensure that the magnitudes of A (l and A s correspond to the 
specified values of P* and P* , respectively, is described in the following section and in 
Appendix A. 

Prescribed edge displacements and rotations 

Edge displacements and rotations are applied to a shell by specifying values for the 
displacements and rotations shown in Fig. 4 that correspond to the positive-valued stress 
resultants shown in Fig. 1. In particular, the axial and tangential displacements, u n and u s , 


respectively, that are applied to the boundary segment, f (/) , are given by 


w 1 (n.ei) = w* 


(8a) 


U 2 [( e 3 Xn ) #e 2 ] = U ] 


(8b) 


Similarly, the transverse displacement u* and the rotation about an axis tangent to an edge 
that are applied to the C b boundary segment are defined by 


w 3 = u t 


(9a) 


n 



M n * e i) = ^ 


(9b) 


As mentioned previously, these prescribed displacements are enforced through the use of 
elastic edge restraints (springs) to relax kinematic admissibility requirments on the functions that 
are used to represent the displacement fields. The uniformly distributed extensional and 
rotational springs that are attached to the shell edges in the normal, tangential, and transverse 
directions and used to enforce the kinematic boundary conditions are depicted in Fig. 4. 

Specifying appropriate stiffness values for the springs results in full or partial restraints along 
the shell edges. A zero value of the spring stiffness corresponds to a traction-free-edge 
condition. In contrast, a value of the spring stiffness that is large compared to the corresponding 
shell stiffness effectively corresponds to a prescribed zero-valued boundary displacement or 
rotation. This approach effectively yields a prescribed kinematic boundary condition in the limit 
as the relative stiffness of the spring becomes much greater than the corresponding shell 
stiffness. Similarly, values for the spring stiffness can be selected that correspond to a specified 
uniform elastic restraint along an edge, similar to that provided by a rigid end-ring. This 
capability is important, and useful, because in some test fixtures or actual structures the edge 
supports may not be stiff enough to simulate a fully clamped boundary condition or flexible 
enough to simulate a simply supported boundary condition. 

As depicted in Fig. 4, the membrane displacements, u n and u s , and the transverse 

displacement, u t = w 3 along the f h boundary segment are restrained by extensional springs with 
stiffness values of S n , S s , and S t in the directions normal, tangent, and transverse to the 
boundary, respectively. In addition to the extensional springs, the edge rotations, 9 n and 9 S , 
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along the I th boundary segment are restrained by rotational springs with stiffness values of J n 
and J s that correspond to rotation about axes tangent and normal to the edge, respectively. 

Extensional springs in the directions nonnal and tangent to the shell edge, with stiffness 
values of s n and.v s , are also used to simulate load introduction through a rigid end-ring, as 

shown in Fig. 3. Specifying values for the spring stiffnesses s n and s s that are relatively large 

compared to the corresponding shell stiffnesses causes the shell edge to behave as if a rigid end- 
ring is attached that produces the uniformly distributed displacements with the corresponding 
magnitudes A n and A s . The values for A n and A. that correspond to the specified concentrated 

loads are detennined by using a penalty parameter approach. This approach enforces the 
difference between the edge displacements of the shell and the unknown unifonn rigid end-ring 
displacements, (w (I -A (i )and (u s -A s ) to vanish, while retaining the corresponding potential 

energy of the applied concentrated loads P* and P* . 

Kinematics and Stress-Strain Relations 

The kinematic equations used in the present study are based, to a large extent, on the 
assumptions of Love-Kirchhoff classical thin-shell theory. Specifically, the axial, 
circumferential (tangential), and nonnal (normal to the mid-surface) displacements of a generic 
point of the shell are denoted byU 1 (s 1 ,s 2 ,s i ) , U 2 (.s ] ,s 2 ,s 2 ) andF,(.s' l ,.s' 2 ,v ( ) , respectively. The 
corresponding displacements of a generic point of the shell mid-surface that share the same unit 
vector normal to the mid-surface are denoted byu l (s l ,s 2 ) , u 2 (s l ,s 2 ) andM 3 (.s l ,.s' 2 ) , respectively. 
In classical shell theory, these displacements are related by 

U l (s l ,s 2 ,s 2 ) — u 1 (s l ,s 2 ) ~ s 3 /t (,V| . ,v. ) (10a) 
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U 2 (s 1 ,s 2 ,s 3 ) = u 2 (s l ,s 2 )-s 3 ^ 2 (s lt s 2 ) 


(10b) 


U 3 (s 1 ,s 2 ,s 3 ) = u 3 (s 1 ,s 2 ) (10c) 

where / 0 1 (s 1 ,s 2 ) and /3 2 (s l ,s 2 ) are the mid-surface rotations about the .s '2 and ,S| axes, 
respectively, that are given by 


ft(s t ,s 2 ) = ii 3J (s n s 2 ) 


(11a) 


p 2 (.s, ,s 2 ) = u 3 2 (.S’, ,s 2 )~ u 2 (,s, , s 2 ) 

R(s 2 ) 


(lib) 


in which a subscript after a comma denotes partial differentiation. The corresponding linear 
membrane-strain-displacement relations are given by 




°22 
y 12 




U n 7 H U-> 

v ’ R v 

( W l,2 U 2,\ ) 


(12a) 


and the bending-strain-displacement relations are given by 


K = « 


1C, 


K 


22 


1C, 


12 


~U 


3,11 

f u 2 ^ A 


vi?y,2 


-2 


1 


(12b) 


It is important to point out that the expression given for the change in surface twist due to 
deformation, K n , is that originally published by Love" ’ in 1888 for general shells, in terms of 
lines of principal-curvature coordinates, and derived in the book by Timoshenko and 
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Woinowsky-Krieger for circular cylindrical shells. As indicated by Bushnell, the expression 
for k X2 vanishes for rigid-body motions in contrast to the corresponding expression presented in 

Reissner's version of Love's first-approximation shell theory (see Reissner, Kraus, and 
Naghdi ). Equations (12a) and (12b), and the more general forms presented by Bushnell,' are 
sometimes referred to as the Love-Timoshenko strain-displacement equations. Justification for 

•50 

this terminology is given by Chaudhuri. 

The stress-strain relations used in the present study are those of the classical theory of 
laminated plates and shells, 34 which are based on a linear through-the-thickness distribution of 
the strain fields. For a thin, symmetrically laminated cylindrical shell, with variable wall 
thickness, the relationship between the membrane and bending stress resultants and the 
membrane and bending strains is expressed conveniently in matrix notation by 


N = A(j 15 j 2 )£ (13a) 

and 

M = D(j 1 ,s 2 )k (13b) 

The membrane and bending stress resultants in Eqs. (13a) and (13b) are defined as 
N r ={V,„V 22 ,V l2 } (14a) 

and 

M T ={M n ,M 22 ,M l2 } (14b) 


It is important to reiterate that when shell-wall thickness variations are present, the membrane 
and bending stiffness matrices, A(s l ,s 2 ) and Diy ,.s\) , are dependent on the curvilinear surface 
coordinates .s’, and.s' 2 . 
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It is convenient, in the present study, to combine the relations given in Eqs. (13a) and (13b) 


into the matrix form 

s = Ce 


in which S, e and C are defined as follows: 
s r ={N r ,M r } 

e r ={ £ r ,K r } 


C = C (s lt s 2 ) 


A(j 19 j 2 ) 0 

0 D(.s',„9 2 ) 


(15) 


(16a) 

(16b) 

(16c) 


Equations Governing the Response 

A general analytical approach for the exact solution of the equilibrium equations for a 
laminated-composite cylindrical shell with variable curvature is not mathematically tractable. 
Therefore, a semi-analytic variational approach that is based on the principle of stationary 
potential energy is used in the present study to obtain numerical results. Because elastic edge 
retraints are used as a means to relax the kinematic admissability conditions on the assumed 
displacement functions, and because a rigid-end-ring capability is used to impose shell-end force 
resultants, the potential energy consists of the elastic strain energy of the shell and the elastic 
edge restraints and the potential energy of the applied loads. In particular, the potential energy is 
expressed symbolically by 

^(q, A) = U (q) + Q (q, A) + F (q, A) (17) 

in which U and Q represent the strain energy of the laminate and the elastic edge supports 
(springs), and V represents the potential energy due to external boundary loads. Their explicit 
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forms are presented in Appendix A. The symbol q is the vector of unknown, generalized 
displacement coefficients that arises from the mathematical representation of the mid-surface 
displacement fields that is used in the variational solution process. In particular, the mid-surface 
displacement fields are given symbolically by Wj(q) , w 2 (q) , andw 3 (q) . The symbol A represents 

the vector of unknown edge displacements that arise from prescribing end loads. 

Subjected to the constraint equations that arise from the use of Lagrange multipliers, the 
equations governing the shell response are obtained by enforcing the requirement that the first 
variation of the total potential energy vanish. As discussed by McFarland et ah, because the 
constraint equations are not functionally dependent on spatial coordinates, y and.v 2 , the 

equations governing the response may be generated by modifying the total potential energy into 
the form 

n (q> A, X.) = ;r(q, A) + W (q, k) (18) 

in which W is viewed as the potential energy arising from constraint reactions. In particular, 

W (q, X) = k T G q = 0 (19) 

where k is the unknown vector of Lagrange multipliers and G is the known constraint 
coefficient matrix. 

Substituting the specific expressions forL(q) , Q(q,A), V (q. A) , and W (q,k) that arise 
from approximation of the surface-displacement field and enforcing the first variation of the 
modified form of the total potential energy to vanish lead to 

Sn = <Jq r [k„ q + S„ q-s, 4 A-f'-T* +G r x] 

+ <SA r [ Sia A-s' a q-P*]+<S>/Gq = 0 (20) 
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in which the matrix, k represents the stiffness matrix of the shell and requires evaluation of the 

corresponding integrand over a doubly connected region (see Appendix A for details). The 
spring-stiffness matrices, S w ands AA , are associated with the deformation of the shell edges and 

displacement of the rigid end-ring, respectively. The spring-stiffness matrix, s ?A , captures the 

coupling between the displacement of the shell edges and the rigid end-ring. The vectors f*,T*, 

and P arise from the prescribed boundary displacements, external tractions and moments, and 
the concentrated forces applied to a rigid end-ring, respectively. For the arbitrary variations 
( 6>'q, SA , and 8 /. ), the stationary condition requires that the following equations must be 
satisfied: 


[K+S t ,)q-s l4 4-f-T-+G r x] = 0 

(21a) 

A-sJ s q-P‘] = 0 

(21b) 

0 

-O 

II 

o 

(21c) 


It is convenient, to express Eqs. (21a) - (21c) into the single matrix equation 

KQ = F 


( 22 ) 


where K and F represent the overall, system stiffness matrix and the overall load vector, 
respectively. These matrices have the general, expanded fonn 



sy 

sy 

a 

1 

1 

-Si 

> 

G T 


F*' 

K = 

-s r 

S AA 

0 r 

and F = < 

p* 


G 

0 

0 


0 


(23a, b) 
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in which 


K «= k « +S « and F* =f* +T* 

The vector of unknowns, Q , that appears in Eq. (22) is defined as 


Q = 


(23c, d) 


(24) 


Solving for the vector of unknowns in Eq. (22) yields all the infonnation needed to obtain a 
complete variational solution to a specific problem. The accuracy of a solution depends on the 
number of terms included in the expressions for the local and global functions representing the 
displacement fields and converges to the corresponding exact solution as the number of terms 
increases. 


Displacement-field representation 

Representation of the mid-surface displacement field is a critical step in the variational 
solution to the problem. By relaxing the requirements for kinematic admissibility, the mid- 
surface displacement fields are represented in the present study by a combination of rigid-body 
modes, u Ri , and global and local functions, denoted by u i and u i , respectively; that is, 

Uj = u Ri + u j +u j (25) 

where the values of the index are given by i= 1,2, and 3. The rigid-body modes account for the 
overall or global translation and rotation of the shell, and are selected so that they produce 
neither membrane strain nor changes in shell curvature and twist. These terms are included for 
the completeness of the kinematics of the cylindrical shell. The presence of the appropriate 
displacement boundary conditions inherently eliminates the rigid-body motion. However, for 
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cases where an insufficient number of kinematic boundary conditions are imposed, these rigid- 
body terms need to be eliminated, as discussed in detail in Appendix C. Following the complex- 
variable solution techniques used in the theory of elasticity, the local functions are expressed in 
tenns of robust, unifonnly convergent Laurent series (used for doubly connected regions) to 
enhance capturing steep stress gradients and deformations near the cutout. Complete sets of 
trigonometric expansions are used to primarily capture the overall global response of the shell. 
Here, completeness means that all the fundamental waveforms needed to construct the typical 
overall deformations of a shell are included in the set. 

For convenience, the displacement representations are rewritten in matrix form as 


U, =V T Ri a R +Vfc i+ V[ a (/ = 1,2) 

(26a) 

U 3 = a R + y[ C 3 + v 3 r p 

(26b) 

An even more useful, compact form is given by 


u i = V/ q with /' = 1,2,3 

(27) 

where the vector of unknown displacement coefficients, q , is defined by 


q r ={a^,c[,C 2 ,c 3 ,« r ,p r } 

(28) 


In Eq. (28), the vector a R contains the unknown coefficients for the rigid-body motion of the 
shell, and the vectors a and p contain the real and imaginary parts of the unknown coefficients 
a„ m and J3 nm , respectively, that are associated with the local functions. The vectors c, , 
where / = 1,2,3 , contain the real-valued unknown coefficients, c i(mn) that are associated with the 
global functions. The explicit forms used herein for the unknown coefficient vectors , C, , a , 
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and P that appear in Eqs. (26a) and (26b) along with the vector functions V, (and the 


corresponding subvectors \ Ri , V, , and V, ) are given in Appendix B. 

In addition to the general representation of the shell surface-displacement fields, similar 
matrix expressions are needed for the displacements and rotations of points on the shell 
boundary. In the present study, the boundary displacement vector u r is introduced that consists 
of the mid-surface boundary displacements in the directions normal, tangent, and transverse to a 
shell edge, and the mid-surface rotations about axes that are nonnal and tangent to a shell edge. 
The boundary displacements in the directions normal, tangent, and transverse to a shell edge are 
denoted herein by a n , u s , andw, , respectively. Similarly, the mid-surface rotations about axes 

that are tangent and nonnal to a shell edge are denoted by & n and l 9 , respectively. In terms of the 

vector of unknowns defined by Eq. (28), the boundary displacements and rotations are expressed 
in matrix form by 

u r =Bq (29) 


in which the boundary displacement vector, u r is defined by 

u r ={ U n’ U s’ U ,A) 

The matrix B is a known matrix of coefficients that is defined as 

u r 

n 

u T 

B = " 

in which the sub-vectors, , u[ , uf and are known and defined by 

u '=( n *e i)Vj 


(30) 


(31) 


(32a) 
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(32b) 


u.v =[(e 3 xn>e 2 ] V 2 

uf=V 3 r (32c) 

and 

9« = ( n * e i) V 3 ^ (32d) 

Strain- and stress-resultant-field representation 

After defining the shell mid-surface displacement field in terms of the generalized coordinate 
q, the corresponding representation of the strains is obtained by substituting Eq. (27) into the 
strain-displacement relations given in vector form by Eqs. (12a) and (12b). This substitution 
yields 

£ = L £ q (33a) 

and 


K = L^q 


where the strain-coefficient matrices L, and L, are defined as 


L = 


\ T 

’u 

\ T + — \ T 
’ 2,2 ' n ’3 

K 

\ T -t-v 7, 
’ 1,2 T ’ 2,1 


-V 

’ 3,11 

I T R „ 

-V + — V +^-V 

* 3,22 R 22 R 2 2 


_2V r +— V r 

^ ’3,12 T r ’2,1 


(33b) 


(34a) 


(34b) 
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Next, the representations for £ and k are substituted into Eq. (15b) to obtain 
e = Lq (35) 

where the overall strain-coefficient matrix L is defined as 

L r =[L ; L r r ] (36) 

Finally, the corresponding matrix representation of the stress resultants in terms of the 
generalized coordinates is obtained by substituting Eq. (35) into constitutive Eq. (15). The 
resulting vector of stress resultants is given by 

s = C Lq (37) 

Constraint Equations 

In the generalized-coordinate representations for u x and« 2 , the coefficients c 1(00 ) and c 2(0 o) 

associated with the global functions, u x and u 2 , also correspond to rigid-body translation in the 

.s', direction and rigid-body rotation about the .s, axis, respectively. These two redundant rigid- 
body modes are eliminated by introducing constraint conditions using Lagrange multipliers. In 
particular, the unknown Lagrange multipliers A KK/j(l) and /. AV(/;(2) are associated with the 

redundant rigid-body modes. Also, multi-valuedness of the normal-direction displacement 
u 3 (s x ,s 2 ) that arises from the presence of logarithmic terms in the Laurent-series-expansion for 

the local function-must be eliminated. The u nkn own Lagrange multipliers A, sv(r) and A sv(s) are 

used herein to eliminate this multi-valuedness. Likewise, the rigid-body modes of the cylindrical 
shell must be eliminated by the Lagrange multipliers h RB( j) ( j = 1, 6 ) if the specified kinematic 

boundary conditions are not sufficient enough to prevent them. In other words, the non- 
vanishing rigid body modes must be eliminated by introducing constraint conditions prior to the 
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stress analysis in order for the overall system stiffness matrix K, given in Eq (22), to be 
nonsingular. 

These requirements on the representation of the shell displacement field are enforced by 
using constraint equations that use Lagrange multipliers. These constraint equations are 
functionally independent, forming a set of linearly independent equations equal in number to the 
total number of Lagrange multipliers. The Lagrange multipliers can be viewed as the reactions 
that are needed to enforce the corresponding constraints. In the present study, all of these 
constraint conditions are included in the matrix equation given in Eq. (19). The explicit form of 
the vector of unknown Lagrange multipliers, k , and the known coefficient matrix, G , are given 
in Appendix C. 

Overview of Validation Studies 

A limited series of validation studies were conducted in the present study to determine the 
accuracy of results obtained by using analysis method presented herein. Specifically, the studies 
included circular and non-circular cylindrical shells with either a circular or an elliptical cutout 
under unifonn tension. The stress resultants around the circular and elliptical cutout for varying 
aspect ratios and orientations in a circular cylinder as well as the stress concentrations arising 
from a circular cutout in a non-circular cylindrical shell were computed. Comparisons of the 
stress-resultant distributions and magnitudes in the shells were made with the corresponding 
results obtained by using an in-house finite element program developed earlier by Madenci and 
Barut. This finite element program has been validated, to a large extent, against previously 
published experimental and numerical results for stress, buckling, and post-buckling of thin-shell 
structures (see Madenci and Barut ’ ). Therefore, this finite element program is expected to 
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serve as a reliable indicator of the accuracy of the analysis methods and results presented herein. 
Overall, the comparisons indicate very good agreement (less than 1% difference) between the 
corresponding results produced by the two analysis methods. For shells with high-aspect-ratio 
cutouts, differences of approximately 5% were obtained and found to be the result of insufficient 
mesh refinement in the finite element models. 

Selected Numerical Results 

Selected numerical results are presented in this section to demonstrate the utility of the 
analysis method presented herein and the potential for its use in developing design technology. 
These results elucidate the effects of loading condition, non-circular cross-section geometry, 
wall-thickness variation, cutout shape, cutout size, and cutout orientation on the intensity of 
stress-resultant concentrations near a cutout. Specifically, tension, torsion, and pure-bending 
loads are considered for [45° /- 45° / 90° / 0° / 90° /- 45° / 45°] s quasi-isotropic shells with length 
L = 356 mm and made of graphite-epoxy plies. The nominal ply thickness isfi. 0 = 0.14 mm , 
resulting in the total thickness of the shell given by h = 2.24 mm , and the ply orientation angles 
are measured with respect to the longitudinal shell axis. The Young’s moduli of each ply in the 
longitudinal, fiber direction and in the direction transverse to the fibers are specified as 
E l = 135. 0 GPa and E T = 13.0 GPa , respectively. The in-plane shear modulus and Poisson’s ratio 

of each ply are given by G LT = 6.4 GPa andv ir = 0.38 . 

The effects of varying the radius of curvature R 0 on the stress-resultant concentration along 
the contour of a circular cutout with radius a = 25.5 mm are shown in Fig. 5 for a circular 
cylindrical shell subjected to a uniform axial tension load. Four curves that correspond to values 
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of R 0 /L= 0.5, 0.75. 1, and 1.25 are presented that show the tangential stress resultant, N H 
nonnalized by the far- Held applied uniform stress resultant A 0 , as a function of position around 
the cutout (indicated by the "cutout angle", (j)). As shown in Fig. 5, the stress-resultant 
concentration is a maximum at ^ = 90° and 270° (at the net section of the shell) for each case 
and reduces from a maximum value of approximately 4.0 to a minimum value of 3.4 at the net 
section as the radius of curvature increases. In addition, the results show that the N^(a, 90°) 

stress-resultant concentration approaches the well-known value of three for an isotropic plate as 
the shell radius increases. Away from the net section, changes in the radius of curvature have a 
relatively small effect on the stress-resultant concentration. 

The effects of varying the circular-cutout radius on the stress-resultant concentration along 
the contour of a circular cutout is shown in Fig. 6 for a circular cylindrical shell with radius 
R 0 =381 mm and subjected to a uniform axial tension load. Five curves that correspond to 
values of the cutout radius a = 15, 25.5. 30, 40, and 50 mm are presented that also show the 
tangential stress resultant N r/ir/) (a,<f>) , normalized by the far- Held applied uniform stress resultant, 

A 0 , as a function of the cutout angle (f) . The results in Fig. 6 show that the stress-resultant 

concentration is a maximum at the net section of the shell for each case, as expected, and 
changes significantly from a minimum value of approxiamtely 3.1 to a maximum value of 5.1 at 
the net section as the cutout radius increases - an increase of approximately 65%. The results 
also show that the N^(a, 90°) stress-resultant concentration approaches the well-known value of 

three for an isotropic plate as the cutout radius decreases. Away from the net section, changes in 
the cutout radius have a much smaller effect on the stress-resultant concentration. 
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The effect of varying the elliptical-cutout aspect ratio, a/b , on the tangential stress-resultant 
distribution around the edge of a cutout in a cylindrical shell with radius R 0 =178 mm, and 
subjected to uniform tension is presented in Fig. 7. The orientation of the elliptical cutout is 
specified by^ = 0° . Two curves that correspond to the locations </> = 0° and 90° are presented 
that show the tangential stress resultant, N„n (a 0 ,j3 ) normalized by the far- field applied uniform 
stress resultant A 0 , as a function of the cutout aspect ratio. As expected, the normalized stress- 
resultant concentration, N^(a 0 ,j3)/N 0 , remains negative for all aspect ratios at (/) = 0" , 

consistent with the expected Poisson effect, and the magnitudes are relatively insignificant at this 
location. In contrast, large stress-resultant concentrations are indicated at the net section (cf> = 
90°) that diminish from a maximum value of approximately 17.0 for a widthwise, slot-like cutout 
with (a = 5 mm and b = 30 mm) or (a/b = 1/6) to a minimum value of 1.4 for a lengthwise, slot- 
like cutout (a = 30 mm and b = 5 mm) or (a/b = 6) . 

The effects of varying the orientation of a high-aspect-ratio, slot-like elliptical cutout on the 
stress-resultant concentration along the cutout contour is shown in Fig. 8 for a circular 
cylindrical shell with radius f? 0 =178 mm and subjected to a uniform axial tension load. The 

major and minor axes of the cutout are given by a = 30 mm and b = 5 mm, respectively. The 
orientation of the elliptical cutout, with respect to the longitudinal shell axis, is measured by the 
angle, i// . Three curves that correspond to values of i// = 0°, 45°, and 90° are presented that show 
the tangential stress resultant at the cutout edge, N „„ normalized by the far-field applied uniform 
stress resultant, N 0 as a function of the cutout angle (/> . 

The results in Fig. 8 show that the stress-resultant concentration is the least pronounced for 
the case of y/ = 0°. For this case, the cutout major axis is aligned lengthwise with the shell axis 
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and the net section of the shell is the largest. The location on the cutout edge defined by <j> = 0° 


corresponds to where the edge of the cutout intersects the major axis. At this location, the edge 
of the cutout is in tangential compression (N^/N 0 =-1.6), consistent with a Poisson effect. 

The location defined by 0 = 90° corresponds to where the edge of the cutout intersects the minor 
axis; that is, at the net section of the shell. At this location, the edge of the cutout is in tangential 
tension (N^/N 0 = 1.4). Between approximately ^ = 10° and 170° and between 0 = 190° and 
350°, the cutout width (and hence net section width) does not vary greatly. This attribute 
accounts for the corresponding flat regions in the y/ = 0° curve shown in Fig. 8. 

For the case of y/ = 90 " , the cutout major axis is perpendicular to the shell axis and the net 
section of the shell is the smallest. As before, the locations defined by </> = 0° and 180° correspond 
to where the edge of the cutout intersects the major axis; that is, at the net section of the shell. 
The results in Fig. 9 show that the edge of this high-aspect-ratio cutout has extremely high stress- 
resultant concentrations at these locations (N /3/3 /N 0 = 17.) that have very step gradients. 

Between approximately 0 = 5° and 175° and between 0 = 185° and 355°, the analysis predicts 
relatively benign variations in the stress-resultant concentration. The case oft// = 45" , exhibits 
stress-resultant concentrations that are, for the most part, bounded by the corresponding results 
for y/ = 0° and 90°. The analysis also predicts very high stress-resultant concentrations where 
the cutout edge intersects the major principal cutout axis ( N pp / N 0 = 8.2 ). 

The effects of varying the cross-section eccentricity (see Eq. (2)) of a tension-loaded oval 
shell with a circular cutout are shown in Fig. 9. The results in this figure correspond to the 
equivalent shell radius R 0 = 381 mm and a circular-cutout radius given by a = 25.5 mm . Moreover, 
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the tangential stress-resultant concentation at the shell net section, tv (a, 90°) , normalized by the 
applied load N 0 , is shown as a function of the eccentricity parameter for the range of 
-0.15 <£<0.15. As indicated in the figure, negative and positive values of £ correspond to 
cylindrical shells with the largest cross-sectional width oriented parallel and perpendicular to the 
tangent plane that passes through the two points of the cutout edge that are on the surface 
generator that passes through the center of the cutout, respectively. A value of £ = 0 
corresponds to a circular cross-section and a value of £ = 0.15 corresponds to cross-sectional 
aspect ratio of 0.9. 

The results presented in Fig. 9 show that the stress-resultant concentration is affected 
benignly by the cross-sectional eccentricity. In particular, the stress-resultant concentration 
increases almost linearly with increases in the eccentricity parameter from N^(a,90°)/N 0 = 3.5 

to 3.6, which is slightly less than a 3% variation. This trend is understood by noting that the 
shells that correspond to negative values of £ are flatter near the cutout than those that 
correspond to positive values of £ and, as indicated by the results in Fig. 5, are expected to have 
the lower values for the stress-resultant concentrations. 

The effects of longitudinal and circumferential periodic wall-thickness variations on the 
stress-resultant concentration at the net section of circular cylindrical shell with 
radius R 0 =178 mm, circular cutout radius a = 25.5 mm, and subjected to uniform axial tension 

load are shown in Fig. 10. Two monotonically increasing curves that correspond to values of s x 
(with s 2 = 0 ) and s 2 (with e x = 0 ) are presented that show the tangential stress 
resultant N j(j (a, 90°), normalized by the far-field applied unifonn stress resultant A n , as a 
function of thickness-variation amplitudes (see Eq.(3)) that range from 0 to 0.2. For the 
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longitudinal thickness variation, the wave numbers used in Eq. (3) are m l =\ and/n 2 = 0 . 
Similarly, for the circumferential thickness variation, the wave numbers used in Eq. (3) are 
m 1 = 0 and m 2 = 1 . 

The results shown in Fig. 10 indicate that the stress-resultant concentration at the shell net 
section increases as the magnitude of the thickness variation increases, for variations in either the 
longitudinal or circumferential direction. The maximum variation in the results is approximately 
56%. Furthermore, the change in the stress-resultant concentration is slightly more pronounced 
for the circumferential thickness variation than for the longitudinal thickness variation. These 
increases are primarily due to a drastic loss of bending stiffness near the net section of the shell, 
as indicated by the wave numbers m 1 = 0 and m 2 = 1 , where the thickness of the shell near the 
center of the cutout is smaller. 

The effects of varying the radius of curvature R 0 on the stress-resultant concentration along 

the contour of a circular cutout with radius a =25.5 mm is shown in Figs. 11 and 12 for a 
circular cylindrical shell subjected to a unifonn torsion load and a pure-bending load, 
respectively. The pure-bending load corresponds to using t n = M 0 cos(;r/s 2 ) in Eq. (5a). Four 
curves that correspond to values of R 0 /L = 0.5, 0.75. 1, and 1.25 are presented that show the 
nonnalized values of the tangential stress resultant as a function of position around the 
cutout. In Fig. 11, N H is normalized by the far- field applied uniform shear stress resultant, T 0 . 
In Fig. 12, A^is normalized by the far- Held applied uniform bending stress resultant, M 0 . 

The results in Fig. 1 1 indicate that the stress-resultant concentration has identical maximum 
magnitudes at (j) = 45°, 135°, 225°, and 315° (at the net section of the shell) for each case, which 
corresponds to maximum diagonal tension and compression stress resultants associated with the 
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shear stress resultants near the cutout. The magnitudes of the stress-resultant concentration for 
these four locations reduces from a maximum value of 6.8 to a minimum value of 5.1 as the 
radius of curvature increases (33% variation). Away from these four locations, changes in the 
radius of curvature have a smaller effect on the stress-resultant concentration. The results in Fig. 
12 indicate that the stress-resultant concentration for the shell subjected to the pure-bending load 
is quite similar to that presented in Fig. 5 for the corresponding tension-loaded shell. 

Specifically, the stress-resultant concentration is a maximum at <j) = 90° and 270° (at the net 
section of the shell) for each case and reduces from a maximum value of 4.0 to a minimum 
value of 3.5 at the net section as the radius of curvature increases (14% variation). In addition, 
Na/M o approaches the well-known value of three for an isotropic plate as the shell radius 
increases, and away from the net section, changes in the radius of curvature have a relatively 
small effect on the stress-resultant concentration. For the case of R 0 /L = 0.5 shown in Fig. 12, a 

contour plot of N Hl /M 0 near the cutout is shown in Fig. 13. The extent of the stress 
concentration at the shell net section (0=90° and 270°) is clearly captured by the analysis 
method presented herein. The highest stress-resultant concentration is N^(a,90°)/M 0 = 4 and 

it attenuates to the value of 1.01 at a radius of about 80 nun (approximately three times the 
cutout radius), measured from the center of the cutout. 

Concluding Remarks 

A special-purpose, semi-analytical approach based on complex potential functions has been 
presented that can be used to investigate the behavior of thin, noncircular cross-section 
cylindrical shells made of laminated-composite materials and with a cutout, efficiently and 
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parametrically. In particular, the effects of radius of curvature; elliptical cutout size, aspect ratio, 
and orientation; oval cross-section eccentrictity; wall-thickness variations; and loading 
conditions on the stress-resultant concentration near the cutout have been presented for a quasi- 
isotropic shell subjected to uniform tension, unifonn torsion and pure bending. In addition, 
studies that were conducted to validate the analysis method have been described. 

A key finding of the results obtained with this analysis method is that the maximum 
tangential stress-resultant concentration near a circular cutout in a tension-loaded, circular, quasi- 
isotropic shell increases by approximately 18% as the shell radius-to-length ratio decreases from 
1.25 to 0.5. Likewise, increases in the maximum tangential stress-resultant concentration as 
large as 65% have been found to occur with a five-fold increase in cutout radius. Results have 
also been presented that show extremely high tangential stress-resultant concentrations can occur 
for high aspect ratio elliptical cutouts whose principal axes are not aligned with the longitudinal 
axis of a tension- loaded shell. 

Additionally, results have been presented that show tension-loaded oval shells with a circular 
cutout on one of the flatter sides exhibit slightly lower tangential stress-resultant concentrations 
than the corresponding shell with the cutout on one of the more highly curved sides. Results have 
also been presented that show that wall-thickness variations in either the longitudinal or 
circumferential directions significantly affect the stress concentration, with respect to that for the 
corresponding shell with a nominal thickness. The analysis also predicts that a quasi-isotropic 
shell with a circular cutout and subjected to pure bending that yields the maximum tensile stress 
resultant at the longitudinal axis of the cutout behaves similarly to the corresponding tension- 
loaded shell. The corresponding shell subjected to torsion was found to exhibit the maximum 
tangential stress-resultant concentrations at locations consistent with the maximum diagonal 
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tension and compression near the cutout. Overall, the results demonstrate that the analysis 
approach is a powerful means for developing design criteria for laminated-composite shells. 

Acknowledgement 

The authors wish to dedicate this paper to the memory of Dr. James H. Starnes, Jr. of the 
NASA Langley Research Center. Dr. Starnes was an internationally recognized expert in 
aerospace structures technology and a proponent of the development of special-purpose, design- 
oriented analysis methods such as that presented herein. 

References 

’Lurie, A. I., “Concentration of Stresses in the Vicinity of an Aperture in the Surface of a 
Circular Cylinder f Prikl. Mat. Mekh., Vol. 10, 1946, pp. 397-406. 

"Lurie, A. I., Statics of Thin-Walled Elastic Shells, State Publishing House of Technical and 
Theoretical Literature, Moscow, 1947. 

-j 

Lekkerkerker, J. G., “Stress Concentration Around Circular Holes in Cylindrical Shells,” 
AIAA Journal, Vol. 10, 1964, pp. 1466-1472. 

4 Van Dyke, P., “Stresses about a Circular Hole in a Cylindrical Shell,” AIAA Journal, Vol. 3, 
1965, pp. 1733-1742. 

5 Ashmarin, I. A., “Stress Concentration Around a Circular Opening in an Orthotropic 
Cylindrical Shell,” Prikladnaya Mekhanika, Vol. 2, 1966, pp. 44-48. 

6 Murthy, M. V. V., Rao, K. P. and Rao, A. K., “On the Stress Problem of Large Elliptical 
Cutouts and Cracks in Circular Cylindrical Shells,” International Journal of Solids and 
Structures, Vol. 10, 1974, pp. 1243-1269. 


33 



7 Guz, A. N., Chernyshenko, I. S. and Shnerenko, K. I., “Stress Concentration Near Openings 

in Composite Shells,” International Applied Mechanics, Vol. 37, 2001, pp. 139-181. 

8 

Van Tooren, M. J. L., Van Stijn, I. P. M. and Beukers, A., “Curvature Effects on the Stress 
Distribution in Sandwich Cylinders with a Circular Cut-out,” Composites: Part A, Vol. 3, 2002, 
pp. 1557-1572. 

9 Tennyson, R.C., “The Effects of Unreinforced Circular Cutouts on the Buckling of Circular 
Cylindrical Shells under Axial Compression,” ASME Journal of Engineering Industry, Vol. 90, 
1968, pp. 541-546. 

10 Starnes, J. H., Jr., “Effect of a Circular Hole on the Buckling of Cylindrical Shells Loaded 

by Axial Compression,” AIAA Journal, Vol. 10, 1972, pp. 1466-1472. 

"Pierce, D. N. and Chou, S. I., “Stresses Around Elliptical Holes in Circular Cylindrical 

Shells,” Experimental Mechanics, Vol. 13, 1973, pp. 487-492. 

12 

'Bull, J. W., “Stresses Around Large Circular Holes in Uniform Circular Cylindrical Shells,” 
Journal of Strain Analysis, Vol. 17, 1982, pp. 9-12. 

IT. 

Zirka, A. I. and Chemopiskii, “Stress Concentration in an Axially Compressed Cylindrical 
Shell of Medium Thickness with an Elliptic Opening,” International Applied Mechanics, Vol. 
10, 2003, pp. 1466-1472. 

14 Liang, C., Hsu C. and Chen W., “Curvature Effect on Stress Concentrations Around Circular 
Hole in Opened Shallow Cylindrical Shell Under External Pressure,” International Journal of 
Pressure Vessels and Piping, Vol. 75, 1998, pp. 749-763. 

"Shnerenko, K. I. and Godzula, V. F., “Stress Distribution in a Composite Cylindrical Shell 
with a Large Circular Opening,” International Applied Mechanics, Vol. 39, 2003, pp. 1323- 
1327. 


34 



16 Hicks, R., “Stress Concentrations Around Holes in Plates and Shells,” Proceedings of the 

Applied Mechanics Conference , New Castle Upon Tyne, 1964, pp. 3-12. 

17 Ebner, H. and Jung, O., 1972, “Stress Concentration Around Holes in Plates and Shells,” 

Contributions to the Theory of Aircraft Structures, Delft University Press, Rotterdam. 

18 . 

Li, Y. W., Elishakoff, I., and Starnes, J.H., Jr., “Axial Buckling of Composite Cylindrical 
Shells with Periodic Thickness Variation,” Computers and Structures, Vol. 56, 1995, pp. 65-74. 

19 Sheinman, I. and Firer, M., “Buckling Analysis of Laminated Cylindrical Shells with 
Arbitrary Noncircular Cross Sections,” AIAA Journal, Vol. 32, 1994, pp. 648-654. 

~ Hyer, M. W. and Wolford, G. F., “Progressive Failure Analysis of Internally Pressurized 
Noncircular Composite Cylinders,” 43 nl AIAA/ASME/ASCE/AHS/ASC Structures, Structural 
Dynamics, and Materials Conference, Denver, Colorado, Paper No. 2002-1403, 2002. 

~ Hyer, M. W. and Wolford, G. F., “Damage Initiation and Progression in Internally 
Pressurized Noncircular Composite Cylinders,” 44 rd AIAA/ASME/ASCE/AHS/ASC Structures, 

Structural Dynamics, and Materials Conference, Norfolk, Virginia, Paper No. 2003-1594, 2003. 

22 

Lekhnitskii, S. G., Anisotropic Plates, Gordon and Breach Science Publishers, Inc., New 
York, 1968. 

" Romano, F. and Kempner, J., “Stress and Displacement Analysis of a Simply Supported 
Non-circular Cylindrical Shells under Lateral Pressure,” PIBAL Report No. 415, Polytechnic 
Institute of Brooklyn, New York, 1958. 

24 Shewchuk, J. R, “Triangle: Engineering a 2D Quality Mesh Generator and Delaunay 
Triangulator,” First Workshop on Applied Computational Geometry, Philadelphia, Pennsylvania, 
pp. 124-133, 1996. 


35 



25 

“ Madenci, E. and Barut, A., “Influence of an Elliptical Cutout on Buckling Response of 
Composite Cylindrical Shells with Non-uniform Wall-Thickness and Non-Circular Cross- 
Section,” 44"' AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials 
Conference, Norfolk, Virginia, Paper No. 2003-1929, 2003. 

26 Love, A. E. H., "The Small Free Vibrations and Deformation of a Thin Elastic Shell," 
Philosophical Transactions of the Royal Society of London, Vol. 179, A, 1888. 

“ Love, A. E. H., A Treatise on the Mathematical Theory of Elasticity, 4th ed., Dover 

Publications, New York, 1944. 

28 

“ Timoshenko, S. and Woinowsky-Krieger, S., Theory of Plates and Shells, 2nd ed., McGraw- 
Hill Book Company, New York, 1959. 

29 

“ Bushnell, D., "Computerized Analysis of Shells - Governing Equations," Computers & 
Structures, Vol. 18, No. 3, 1984, pp. 471-536. 

~ Reissner, E., "A New Derivation of the Equations for thr Deformation of elastic Shells," 
American Journal of Mathematics, Vol. 63, 1941, pp. 177-184. 

3 1 

Kraus, H., Thin Elastic Shells - An Introduction to the Theoretical Foundations and the 
Analysis of Their Static and Dynamic Behavior, John Wiley and Sons, Inc., 1967. 

‘Naghdi, P. M., "Foundations of Elastic Shell Theory," Office of Naval Research, Technical 
Report No. 15, January 1962. 

Chaudhuri, R. A., Balaraman, K., and Kunukkasseril, V. X., "Arbitrarily Laminated, 
Anisotropic Cylindrical Shell Under Internal Pressure," AIAA Journal, Vol. 24, No. 11, 1986, pp. 
1851-1858. 

34 Jones, R. M., Mechanics of Composite Materials, 2nd ed., Taylor & Francis, Inc. 
Philadelphia, Pennsylvania, 1999. 


36 



3:i McFarland, D., Bert L. Smith, B. L. and Walter D. Bernhart, W. D., Analysis of Plates, 
Spartan Books, New York, 1972. 

36 

Madenci, E. and Barut, A., “A Free-Formulation Based Flat Shell Element for Non-Linear 
Analysis of Thin Composite Structures,” International Journal for Numerical Methods in 
Engineering, Vol. 37, 1994, pp. 3825-3842. 

Madenci, E. and Barut, A., “Pre- and Postbuckling Response of Curved, Thin Composite 
Panels with Cutouts under Compression,” International Journal for Numerical Methods in 
Engineering, Vol. 37, 1994, pp. 1499-1510. 

Madenci, E. and Barut, A., “Thermal Postbuckling Analysis of Cylindrically Curved 
Composite Laminates with a Hole,” International Journal for Numerical Methods in 
Engineering, Vol. 37, 1994, pp. 2073-2091. 

Appendix A 
Strain Energy of shell 

Based on classical laminated shell theory, the strain energy of the shell can be expressed as 
t/ = |js r e<£4 (38) 

in which A is the planform area of the shell mid-surface. Substituting the expressions for the 
resultant stress and strains, given in tenns of the vector of unknown displacement coefficients, 
q , by Eqs. (35) and (37), leads to 

t/(q) = Ijq r (L r C L) qdA (39) 
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The matrix L involves the derivatives of the assumed, functional displacement representations, 
and C is the overall constitutive matrix defined by Eq. (16c). The expression for the strain 
energy is rewritten into the final form used herein as 

^(q)=^q r k„q (40) 

where 

k () =J(L r CL)<M (41) 

A 

The evaluation of this area integral is performed numerically by employing basic quadrature 
techniques. In this analysis, the quadrature points are pre-determined by employing standard 

24 

triangulation of the entire domain as described by Shewchuk. 

Strain energy of elastic restraints 

The strain energy of the elastic edge restraints (springs), Q , is expressed as 
= £ 1 S.(u.-ufdT 

^ l=\ a=n,s,t r (<) 

S •'„(%- d r ( 42 ) 

a=n,s r (0 

+ Z J S A U a- A a ) 2 dr + 

1= 1 a =n,s r (() 

As depicted in Fig. 4, the boundary displacements u n ,u s , and u t along the l th boundary segment 
are restrained by extensional springs with the stiffness values S n , S s , and S t , respectively. 
Likewise, the boundary rotations <9 n and 3 S are restrained by rotational springs with the stiffness 
values J n and J s , respectively. 
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In order to apply concentrated forces along the edge of a shell and introduce edge 
displacements that are similar to those introduced by a rigid end-ring or by the loading platens of 
a testing machine, additional springs are uses to simulate the load-introduction effects of a rigid 
end-ring. In particular, rigid-end-ring loads are introduced into the shell by using extensional 
springs in the directions normal and tangent to the boundary with corresponding stiffness values 
ofs ; , and.s; v , as shown in Fig. 3. By specifying relatively large values for the spring stiffnesses 

s n and.v s , the laminate edge behaves as if a rigid end-ring is attached that produces the uniform 
displacements A n and A v . In contrast, a relatively small spring stiffness between the shell edge 

and the rigid end-ring eliminates the presence of a rigid end-ring. 

The desired form of the elastic-restraint strain energy is obtained in terms of the unknown 
vector q by substituting expressions for the boundary displacements and rotations, given 
collectively by Eq. (29), into Eq. (42). This step yields 

ndtl (q r s™q+n«;-2q r fr) 

^ 1 = 1 a=n,s,t 

T) 

^ i=\ a=n,s 

f \ 


All 

^ i=l a=n,s 


q r sL«q+ j s a A;</r-2q r s^A c 


-(<) 


where the matrices and represent the stiffness contribution of the extensional and 
rotational springs attached to the i th segment of the boundary. These matrices are defined as 
s 12 = j u a u r a dF (a = n, s, t ) (44a) 


r (0 


and 

J i!=f ( a = n,s ) (44b) 

r w 


(43) 
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The matrix , representing the stiffness of the springs attached to the rigid end-ring, is defined 
as 

S 1X JX u « u X r (a=n,s) (45) 

The load vectors, f^' 1 and , are associated with the prescribed boundary displacements and 
rotations and are defined as 

= j S a u* U„ dr ( a = /1, 5,0 (46a) 

r (0 

and 

■T = | J. K 0„ dr (a = n,s) (46b) 

r (0 

The vector, s^°, is associated with the unknown end-displacements that correspond to a given 
concentrated load and is defined as 

s a )= |X u X r ( a = n,s ) (47) 

4 

The strain energies in the springs that arises from the known prescribed displacements 
( a * , u* and u * ) and rotations ( 9* and 9* ) are defined as 

^u(a) = j S a u a dT (a = n,s,t) (48a) 

r m 

and 

n< £a) = j r ( a = n,s ) (48b) 

r (0 

For convenience, the expression for the strain energy in the springs is recast in matrix form as 
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( 49 ) 


1 


1 


Q (q>A) = -q S qq q + -A s AA A 
-q r V A-q r f* +Q* 

in which the matrices, S , S AA and s A represent the stiffness of the springs associated with the 


deformation of the laminate, the end-displacements and their coupling, respectively. These 
matrices are defined by 


s w = X X s i«+XX J -+XX S 


w 

aa 


t—\ a=n,s,t 


Saa =Diag[^ !) , 


t—\ a=n,s 


«(2) „(1) „(2) 
^ n 9 9 ^ s 


i=l a=n,s 


x2 7i R 0 


e (2) e® e (2 ) 

9 qA 




(50a) 

(50b) 

(50c) 


The vector of unknown end-displacements, A , is defined by 

a t = (a (1) , a <2) , a (1) , a (2) | 

( n 9 n ? s J 


(51) 


The load vectors arising from all prescribed boundary displacements and rotations, f * , is defined 


as 


* 1 \ 1 1 \ 1 *(0 

f = tX X + TXX r ~ 


(52) 


l=\ a=n,s,t 


t= 1 a=n,s 


and the strain energy of all the springs due to prescribed displacements and rotations is 




^ l=\ a=n,& 


7 ^ *— 

^ (,=\ a=n 


Q 


(O * 
9 (a) 


(53) 


Potential of external loads 
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The potential energy of the external tractions (t* n ,t s and t*) and moments (m* n and m* s ) acting 
along the £ th boundary segment, and the concentrated loads (P* and P*) acting on the rigid end 
rings, is given in terms of the corresponding boundary displacements and rotations by 


v =-2 Z I 

l=\ a=n,s,t r (<) 

-ZZ 1 


(54) 


i-\ a=n,s p 


£=\ a=n,s 


Substituting the expressions for the boundary displacements and rotations, given in terms of 
the vector q, and combining terms in Eq. (55) yields 

F(q,A) = -q T’-A P* (55) 

where the vector A, containing the uniform end-displacements A n and A s of the t h boundary 


segment, is defined by 

a t = Ia (1) , a <2) , a (1) , a (2) | 

The load vectors, T*andP* are defined by 

T '' = Z Z I <»+tl | 


(56) 


=1 a=n,s,t p 


t—\ a=n,s p 


and 


p* r _fp*( 1 ) D*( 2 ) p*(l) p*(2) 
r \ r n ’ r n ’ r s ’ r s 


(57b) 


in which P* , with (a = n,s ) , represents the membrane forces applied on the £ boundary 
segment through a rigid end-ring. 


Appendix B 


(57a) 
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Rigid-body modes 

As given by Madenci and Barut 24 , the rigid-body displacements (u m ,u R2 and u R3 ) of a 
cylindrical shell, defined with respect to the curvilinear coordinates, (s x ,s 2 ,s 3 ), are expressed 
herein as 


= a, - a 6 y + a s z 

(58a) 

= a 2 cos 6 - a 3 sin 6 - a 4 ( y sin 0 + z cos 9 ) 
+ a 5 x sin 9 + a b x cos 9 

(58b) 

= a 2 sin 9 + a 3 cos 9 + a 4 (y cos 9 - z sin 9 ) 
- a 5 x cos 9 + a 6 x sin 9 

(58c) 


where 9 denotes the angle between the radius of curvature at a point on the shell surface and z- 
axis as shown in Fig. 1. 


Global functions 

The global functions u i that are used to capture the overall deformations away from the 
cutout are expressed in terms of a series expansion of orthogonal functions of the form 

M m 

U l (S X ,S 2 ) = Y J Yj C iCnn) T m( S l) W n( S 2 ) ( 59 ) 

m= 0 n = 0 

The symbols c i(mn) are the unknown real-valued coefficients, and T m (.s', ) and W n (,s ' 2 ) are defined 
as 


T m (s ,) = 




sin 


(m - 1) 


(C + l) 


m = 0 
m = 1 

m > 1 


(60a) 
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and 


K(* 2 m 


Jcos(«#/2) n=0,2,4,6,8,--- 

jsin((w + 1)0/2) n=l,3, 5,7,9, ••• 


(60b) 


in which -1 < £ <1 and .v, is related to C as .v, = tfL/2, with L being the length of the cylinder. 


Note that W n is periodical. These particular functions were chosen because they fonn a complete 

set of functions when used with Eq. (59). Hence, they are desirable for employing in energy 
based semi-analytic solution techniques such as the total potential energy principal that is used in 
this study. 


Local functions 

The local functions are expressed in terms of mapping functions that transform the contour of 
an elliptical cutout to a unit circle. These mapping functions are used permit the use of Laurent 
series expansions as local functions, which is desirable because Laurent series are analytic and 
uniformly convergent in domains with a circular hole. As a result, the use of mapping functions 
reduces the number of terms in the Laurent series significantly that are needed to adequately 
capture steep stress and strain gradients and local deformations near a cutout. In accordance with 
the principle of minimum potential energy, the local local functions are not required to satisfy the 
traction boundary conditions at the cutout boundary. Thus, the local functions, u i , are expressed 
in the form of Laurent series, in terms of complex functions, as 


Mj =2 Re 


2 N 

y « (i) y a o* (z ) 

m nm nm v sm / 

m = 1 n=—N 

n* 0 


Hip) 


(61a) 
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(61b) 


u 2 = 2 Re 


V u {2) Y « (E>* (z ) 

m nm nm \ em / 

m=l n=—N 
n* 0 


#(/>) 


i/ 3 =2 Re 


2 AT 



m=l n=—N 
n* 0 


#(p) 


(61c) 


with 

P = V x i 2+x 2 (62) 


where the parameter /V defines the extent of the complex series. In these series, and 
are the unknown complex coefficients that appear in Eqs. (26)-(28). The auxiliary function 
H(p ) that defines the domain of influence of the local functions is expressed in a polynomial 
form as 


H{p) = 


1-10 

0 


r^v 

\Po J 


+ 15 


\PoJ 


-6 


ivf 

\PoJ 


0<P<P o 
P> Pc 


with 


H(p o ) = H\p o ) = H\p o ) = 0 


(63b) 


where the prime marks denotes differentiation with respect to the variable p and the parameter 
p 0 denotes the radius of the region in which the local functions are effective. The purpose of 
chosing the auxiliary function is to prevent any possible linear dependency between the local and 
global functions and to restrict the influence of the local functions to a limited domain around the 
cutout. 

The complex functions u ( ^{z £m ) and u (2) (z £m ) that appear in Eqs. (61a) and (61b) are 
defined as 


(63a) 
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ul\ z sJ = ™W Pm( Z sm)~ S c lJ z ,J (64a) 

^\ z em) = ^ys p m (z £m ) + cosy/ q m ( z £m ) (64b) 

where the complex constants p m and q m are given by 

Pm = a x xt*L + a \2 ~ a X(,Psm ( 65a ) 

q m = a n Mem + a 22 1 Mem ~ Q 26 ( 65b ) 


In Eqs. (65a) and (65b), the unknown complex constants, ju em , are the roots to the characteristic 
equation associated with membrane deformation, i.e., 

XPL ~ 2a X6PL + ( 2 «26 + a 66 )PL .... 

. ( 66) 

- 2a 2 6 Psm +a 22 =° 

in which the coefficients a tj are the coefficients of the flexibility matrix a , which is the inverse 
of the stiffness matrix A defined by Eq. (13a). Both the flexibility and the stiffness matrices, 
a and A, are measured with respect to the local coordinate system (jq , x 2 ) . The angle, 
y/ represents the orientation of the local coordinate system with respect to the global coordinate 
system, (sj,s 2 ). 

The complex potential function, ® * nm (z £m ) , appearing in Eqs. (61a) and (61b) is defined as 

<»(0 = C (67) 


in which the mapping functions, ^ £m , map a cutout onto a unit circle. The mapping functions for 

22 

an elliptical cutout, introduced by Lckhnitskif , are given by 



(m = 1 , 2 ) 


( 68 ) 
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where z em = x x + ju sm x 2 , a and b are the major and minor axes of the elliptical cutout, 


and/ = V-l . The sign of the square-root tenn is chosen so that|c,.„, | > 1 (i.e., the mapped point is 

guaranteed to be on or outside the unit circle). 

Inverting the mapping function provides co £m (^ £m ) as 

z =a> (<E ) = r <E (69) 

sm sm^^sm/ sm~sm «> \ / 

^ sm 

in which 


r ~ = \{ a -nbJ ) » S sm=\{ a +i Ms,„ b ) 


(70a,b) 


The unknown complex constants ju sl and /u s2 , and their complex conjugates, i.e., 
ju 6 , 3 = j2 sl and ju eA = /u s2 , are the roots obtained from the characteristic equation associated with 
membrane deformation 

The complex potential functions, F* m (z Km ) in E. (61c) are defined as 


Fnm ( ^ Km ) 


’ Km p. n 0 Km gn — 2 I I . 1 

T 1 /cm 1 5 Km 9 \ n \ ^ 1 

n + 1 n - 1 


-y n = 1 


(71) 


';>L+yC. « = -> 


in which the expressions for the mapping function t; Km and the constants r Km and.sy m have the 
same form as the corresponding expressions for ^ sm , r £m , and s em given by Eqs. (68) - (70a, b), 
except that the subscript s is replaced by k . 

The complex variables z Km are defined by 

**m=*l + /W*2 ( 72 ) 
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in which the unknown complex constants ju kX and ju k2 and their conjugates, i.e., 


fi K 3 = Ji Kl and // a4 = ju h 2 , are the roots obtained from the characteristic equation associated with 
the bending equilibrium equation 

D 22 mL + 4D 2 (,K,„ + (2D n + 4D 66 )j/ Km 

+ ^ D l6^m +D n =0 

where D (/ are the components of the bending stiffness matrix D (see Eq. (13b), which is defined 

with respect to the local coordinate system, (jq ,x 2 ) . 

It is important to note that the local functions in Eq. (61) satisfy the in-plane and bending 
equilibrium equations of a homogeneous, flat laminate of uniform thickness, not a cylindrical 
shell. Therefore, the roots to the characteristic equations, Eqs. (66) and (73) serve as 
approximation to their exact values which are not mathematically tractable. Because the solution 
procedure is based on the principle of minimum potential, their exact values are not necessarily 
required. However, they capture the stress concentration and local deformation near the cutout 
in cylindrical shells because these functions possess the inherent solution characteristics. They 
satify the equilibrium equations exactly as the radius of curvature approaches inifinity and are 
uniformly convergent in a doubly connected region. 

In the displacement representations defined by Eqs. (26a) and (26b), the vectors, a R , c,. , a , 


and P are defined as 


®R 


® R2 ’ ^R3 ’ ^R4 ’ ^R5 


’ a R(>\ 


J _ ( 

' i ~ | C ;(00)’ 


C /(10)’ C i(01)’""’ C !'(M0) 


5 1)1) >••••> C ;(1(M-1)) 5 C i(0M) J 


T | T T T T T T 

® _ |®-V 5®-JV+l’ ,®_1,®1 5 


(74a) 

(74b) 

(74c) 
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in which 


<=K„«y < 75 > 

with 

K = { R e \a nj ] , Im [a nJ ]} (76a) 

and 

in which 

Pi=iPi„py (77) 

with 

Pt={ R e[y]. Im[Aj]} (78) 

The vector functions, V, , with i = 1,2,3 associated with the unknown generalized coordinates q 
appearing in Eq. (27) are defined as 

Vr={vJ,V 1 r ,F,O r ,^,F} (79a) 

V 2 r ={v s r ,,0 r ,V, r ,0 r ,V/,r| (79b) 

v 3 r = {v£, ,0 r ,0 r , Vj r ,! r , Vj r j (79c) 

where 

0 r ={0,0,0,..,0} of order [(M+l)(M+2)/2] (80a) 

0 r ={0,0,..,0} of order 8N (80b) 

The vectors associated with rigid body motion are 
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Vj, ={1,0, 0,0 ,z,-y] 


(81a) 



0, 


dy dz 
ds 2 ’ ds 2 ’ 


f 

y 

V 


dz dv N 

z^— 

ds 2 ds 2 y 


dz 





dz dy 
ds 2 ds 2 


dv dz ^ 

y + z 

ds 2 ds 2 




Similarly, the vectors associated with the global functions are 

\, T ={T 0 ( Sl W 0 (s 2 )JMW 0 (* 2 l 
tmwMjmwm,.. 

T x (jj )W M _ X (s 2 ), T 0 (jj )W M (j 2 )} 


(81b) 


(81c) 


(82) 


in which the expressions for T^) and W t (s 2 ) are given by Eq. (60), and 


\ T = { V r V 

T i ) ’ il-N )’ ’/ 


V 

i(-iV)’ T i(-JV +1) ’i(-l) 

v r y 7, y T 

? v /(!)?•••? T j(7V -1)’ J i(N ) 


(83) 


with 


(84) 


in which 


i{nj) 


2 Re 




,-2Im 




(85) 


where i, j = 1,2 and 
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( 86 ) 


V/= V: 


'3 (-N ) 


T y T yT 

’ ’3 (~N +1)’""’ V- 


V r V 


3(-l) 

T 

3 (TV -1)’ t 3(AT ) 


with 



(87) 


in which 


^-,=!2Re 


F" 


, -2 Tm 


77 * 


(7=1,2) 


( 88 ) 


Appendix C 


Constraint Equations 

The unknown vector of Lagrange multipliers, k and the known coefficient matrix, G , in Eq. 


(19) are defined by 


^ ’ ^RRB( 2) ’ ^SV(r) > ^SV(s) ’ ^RB( 1) ’ \ 


RB( 6) 


and 


o* 

SilfiLB(l) 

o r 

¥ 


(F 

o* 

o r 

%>RRB(2) 

¥ 


IF 

o* 

o r 

¥ 

¥ 

IF 

§ SFM 

o* 

01 

¥ 

o r 

It 

g SV(r) 

855(1) 

o r 

¥ 

o r 

F 

¥ 

T 

&RB( 2) 

o r 

¥ 

o r 


¥ 

T 

855( 3) 

o r 

¥ 

¥ 

¥ 

¥ 

8/25(4) 

o r 

¥ 

¥ 

IF 

¥ 

8/?5(5) 

o r 

¥ 

¥ 

IF 

¥ 

8/25(6) 

o r 

¥ 

¥ 

If 

¥ 


(89a) 


(89b) 
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in which the constant coefficient vectors, g /UM(1) and g^^ are associated with redundant rigid- 

body modes, g SV{r) and g V(/(s) with single valuedness of the radial displacement component, and 

g RB(j) with the rigid-body modes introduced by the global functions defined in Eq. (59). The 

rigid-body modes must be eliminated in the absence of a sufficient number of specified 
kinematic boundary conditions. These terms, as well as, the vectors with zeros are defined in the 
following subsections. 


Redundant rigid-body modes 

The coefficients, c 1(00) and c 2(00) in Eq. (59) for the global displacement functions produce 
additional rigid-body translation in the s l direction and rigid body rotation about the s i axis, 
respectively. Because these rigid-body modes are already represented by a Rl and a R( in Eq. 
(58), the redundant rigid-body motion arising from the presence ofc 1(00) and c 2(00) , must be 


eliminated in order to obtain a unique representation of the displacements. These redundant 
rigid-body modes are eliminated by using the constraint conditions 


^RRB( 1) C 1(00) 0 


(90a) 


^RRB( 2) C 2(00) ^ 


(90b) 


in which the unknown Lagrange multipliers are denoted by a RKH{ , } and h Rm2) . In terms of the 


vector of unknowns, q , these constraints are rewritten in vector fonn as 
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a R 

Cl 

c, 

= 0 

C 3 

a 

P 

C 1 

2 Uo 

C 3 

a 

P 

in which the vectors of zeros, 0 2 , 0 2 and 0^ are defined as 
0* ={0,0,0,0,0,0} (92a) 

0 r ={0,0,0,.., 0} of order [(M+l)(M+2)/2] 

5 r ={0,0,0, ..,0} of order 8N (92c) 

The constant coefficient vector, g RRB is defined as 
g T RRB ={1,0,0,..,0} of order [(M+l)(M+2)/2] 


A, 


'RB( 2) 


0 r 0 r g r o r o r o 7 

^ R-> u ? ORRB ’ V U 


A, 


'RB(l) 


0 T (\T f\T rxT rxT 

R ’ &RRB "> ® ^ 5 0 i 0 


(91a) 


(91b) 


(92b) 


(93) 


Single-valuedness of the radial displacement component 

The multi-valuedness of the normal displacement component that arises from the presence of 
logarithmic tenns in the local expression for the radial displacement component in Eq. (61c) 
must be rendered single valued in order to obtain a unique solution. The logarithmic terms 
associated with complex constants, (3 ]m and f3 hn , with (m = 1,2), in Eq. (71) result in two real 
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constants. Representing the complex variable of the Laurent series, = p m e l 
valuedness requirement is enforced as 

u (4 = pJ e,n ) - (4 = p m e , ' ( ^ +2 " ) ) = 0 

Associated with the complex constants, /3 lm and /? lm with (m = 1,2) , in Eq. (71), 
yields 

m=l 

or 

Im Zk ^ij = 0 

m=l 

and 

2RefEK«. =py-'l 

m = 1 


or 


lm Z[ v ™ A»] =0 


In order to ensure single-valuedness, these constraints are enforced as 


^SV(r) 


y fr B 

L Km . 


= 0 


_m = 1 


and 


^SV(s) 


[“Am Am. 


= 0 


m = 1 


” , the single - 
(94) 

this condition 
(95a) 

(95b) 

(96a) 

(96b) 

(97a) 

(97b) 
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in which the unknown Lagrange multipliers are denoted by A SV(r) and A SV(s) . 

In terms of the vector of unknowns, q, these constraint conditions can be recast in matrix 
form as 




a. 


= 0 


a 




a, 


= 0 


a 

P 


(98a) 


(98b) 


where the constant coefficient vectors g SV(r) and g SF(s) are given by 

T ( T T T 

SSVU) ~ \8sFO)(-A r )’8sKO')(-V+l)’”"’§5K(7)(-l) 

T T T \ 

’ Ssv(j)(l) 1) ’ S^KO'XiV) J 


in which 

gUxi) ={ Im k]» Re [ / i]» Im h]» Re h]} ( 100a ) 

g sv(rXn) = {0, 0, 0, 0} if n 1 (100b) 

and 

g^oxD = { Im [ s i ] > Re [ s i ] > Im [^2 ] , R e [^2 ]} ( 1 00c) 

gsr ( s X n) = ^°’ °’ °’ 0} if n * l ( 100d ) 

with-/V <n< N . 


(99) 
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Rigid-body modes 

In the absence of kinematic boundary conditions, the rigid-body modes of the displacement 
field are eliminated by enforcing the constraint conditions in the form 

^RB(j) a RU) = 0 ( 101 ) 


where j = 1,2, 3,4, 5,6 , and the unknown Lagrange multipliers are denoted by /t AV;( ;) . In terms of 


the vector of unknowns, q , these constraints are rewritten in vector fonn as 


4,0-) gLa).o r ,o r ,o r ,o r ,o I 


a. 


a 


= 0 


( 102 ) 


where the constant coefficient vectors, g RB(j) are defined as 

S r RBU) = { S lj’ S 2jAj’ S *jAj’ S 6j} ( 103 ) 

in which 8 tj is the Kronecker delta. 
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Fig. 1 Geometry, coordinates systems, and applied edge tractions for an oval cylindrical shell 
with an elliptical cutout and non-uniform wall thickness. 


L 

S 5 



Fig. 2 Computational domain of a cylindrical shell with an elliptical cutout. 
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Fig. 3 Application of a concentrated force through a rigid end-ring by using elastic springs 
along a shell edge 



Fig. 4 Types of unifonnly distributed elastic spring supports that can be prescribed along the 
edge of a cylindrical shell. 
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Fig 5 The effect of varying shell radius on the tangential stress resultants around a circular 
cutout in a quasi-isotropic circular cylindrical shell subjected to a uniform tension load. 



Fig 6 The effect of varying circular-cutout radius on the tangential stress resultants around a 
circular cutout in a quasi-isotropic circular cylindrical shell subjected to a uniform tension 
load. 
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1 2 3 4 5 

Elliptical cutout aspect ratio, alb 


Fig 7 The effect of varying elliptical-cutout aspect ratio on the tangential stress resultants 
around a cutout in a quasi-isotropic cylindrical shell subjected to a uniform tension load. 



Fig 8 The effect of varying elliptical-cutout orientation on the tangential stresses around the 
cutout in a quasi-isotropic cylindrical shell subjected to a uniform tension load. 
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Fig 9 The effect of shell cross-section eccentricity on the stress-resultant concentration in an 
oval quasiisotropic cylindrical shell with a circular cutout and subjected to a uniform 
tension load. 



Fig 10 The effects of longitudinal ( e x * 0 and s 2 = 0 ) and circumferential ( £•, * 0 and s 2 = 0 ) wall 
thickness variations on the tangential stress-resultant concentration around a circular 
cutout in a quasi-isotropic circular cylindrical shell subjected to a uniform tension load. 
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Fig 11 The effect of varying shell radius on the tangential stresses around a circular cutout in a 
quasi-isotropic circular cylindrical shell subjected to a unifonn torsion load. 



Fig 12 The effect of varying shell radius on the tangential stresses around a circular cutout in a 
quasi-isotropic circular cylindrical shell subjected to a pure-bending load. 
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Fig 13 Stress resultant distribution near the cutout in a quasi-isotropic circular cylindrical shell 
subjected to a pure-bending load. 
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